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Section  I 


Introduction 


The  Structural  Vibration  and  Acoustics  Branch  of  the  US  Air  Force  Flight  Dynamics 
Laboratory  (YVRDC/FIBG)  has  initiated  a  program  to  study  the  dynamics  and  control  of 
Large  Space  Structures  (LSS).  This  Large  Space  Structures  Technology  Program  (LSSTP) 
is  intended  to  enable  the  Flight  Dynamics  Laboratory  to  instrument,  test,  and  analyze 
large  space  structures  on  the  ground  in  order  to  predict  their  behavior  in  space.  Since  the 
testing  i°  to  be  don^  in  a  ground  based  laboratory,  i.e.  under  1-g  acceleration,  the 
experiments  must  be  designed  so  as  to  counteract  this  gravitational  effect.  One  proposal  is 
to  use  soft  suspension  systems.  Long  cables  can  provide  pendulum  support  with  low 
frequency  for  horizontal  motion,  but,  for  a  pendulum,  the  restraint  in  the  vertical  direction 
is  rigid.  One  way  to  provide  soft  restraint  vertically,  while  maintaining  the  pendulum 
approach,  is  to  counter-balance  the  test  model. 

The  simplest  counter-balanced  suspension  system  is  the  classic  Atwood's  machine.1 
This  consists  of  two  masses  m^  and  m2  connected  by  an  inextensible,  flexible  wire  of 

negligible  mass,  draped  over  a  frictionless,  massless  pulley,  which  in  turn  is  rigidly 
suspended  from  an  overhead  support.  The  motion  of  the  two  masses  is  assumed  to  be 
constrained  to  the  vertical  direction:  either  there  is  no  motion,  the  situation  that  occurs 
when  the  two  masses  are  equal  and  there  is  no  initial  velocity,  or,  as  one  mass  rises  the 
other  falls.  The  most  ambitious  counter-balanced  system  envisaged  in  this  report  consists 
of  a  lattice-type  flexible  structure  suspended  over  pulleys,  no  longer  considered  frictionless, 
by  a  set  of  wires,  no  longer  inextensible,  and  counter-balanced  by  an  identical  structure. 
Since  the  two  structures  do  counter-balance  each  other,  the  net  gravitational  effect,  at 
least  at  the  points  at  which  the  wires  attach  to  the  bodies,  is  null. 

As  might  be  imagined,  there  is  a  host  of  possibilities  between  these  extremes.  In  this 
report  we  record  some  of  these  and  indicate  how  the  motion  of  each  may  be  analyzed.  We 
have  not  studied  in  depth  all  systems  to  be  described;  for  those  that  we  have,  greater  detail 
will  be  given. 

The  theories  and  methods  used  come  from  Lagrangian  mechanics,  linear  elastic 
theory,  the  calculus  of  variations,  non-linear  differential  equations,  numerical  methods  for 
the  solution  of  algebraic-differential  equations,  and  perturbation  theory.  Where  references 
are  appropriate,  these  will  be  offered;  the  report  otherwise  reflects  the  thoughts  of  the 
author. 
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Section  II 


Statement  of  the  Problem 


II.  1. a  The  classical  Atwood's  machine  is  depicted  in  Figure  l.a.  The  analysis  of  this 
device  is  included  only  because  it  serves  as  a  guide  to  the  logic  and  methods  used  in 
subsequent  problems. 

Using  the  symbols  of  Figure  l.a,  the  kinetic  energy  is: 


T  =  (i)m1z12  +  (I)m2z^ 


and  the  potential  energy  is: 


V  -  m1gz1  +  m2gz2 


(1) 


(2) 


The  Lagrangian  function  is: 


L  =  T  -  V  (3) 

There  is,  however,  a  constraint,  namely,  that  the  length  of  the  wire  joining  the  two  masses 
is  constant: 


$  =  (h-Zj)  +  7rR  +  (h-z2)  =  1°  =  length  of  wire  (4) 

The  modified  Lagrangian  function  is: 

L  =  T  -  V  +  A(t)$  (5) 

and  where  A(t)  is  a  function  of  the  time,  t. 

Hamilton's  principle  implies: 

/l 

S  f  L  dt  =  0  (6) 

*0 

which  in  turn  yields  the  Euler-Lagrange  equations: 

m^Zj  +  m^g  -  A  =  0  (7) 

m9z,7  +  m9g  -  A  =  0  ($) 

to  which  we  must  adjoin  the  constraint  condition:  $  =  f . 


■) 


With  this  done,  we  eliminate  A  between  Equations  (7)  and  (8)  and  obtain,  using  Equation 

(A): 


(9) 

(10) 


When  these  equations  are  integrated,  we  obtain: 


2iW  =  ! 


gt2  +  z1(0H  +  Zj(0) 


gt2  +  z2(0)t  +  z2(0) 


(ID 


(12 


■» 


The  values  for  z^(0)  and  z.?(Q)  are  quite  arbitrary  and  depend  on  the  configuration,  i.e.  the 

height  h,  the  radius  R  and  the  length  of  the  connecting  wire.  The  values  for  z^(0)  and 

z9(0)  are  not  arbitrary  since  z^(0)  +  z9(0)  =  0,  i.e.  as  one  mass  rises  the  other  must  fall  at 
the  same  speed,  both  initially  and  for  all  future  time. 


Il.l.b  The  modified  Atwood's  machine  described  in  this  section  is  depicted  in  Figure 
l.b.  The  method  of  solution  parallels  that  given  above  for  the  classical  case. 

The  kinetic  energy  is: 


and  the  potential  energy  is: 


T  =  (^nVi2+  (^)m2z22 

(1) 

V  =  mlgzl  +  m2gz2 

(-) 

We  could  proceed  and  introduce  a  constraint  equation;  it  would  be: 

ZjW  -  z^O)  z9(t)  -  z9(0) 


IT 


+ 


=  0 


(3) 


For  variety,  we  circumvent  the  use  of  the  constraint  by  reducing  the  expressions  for  the 
kinetic  and  potential  energies  (written  above  as  functions  of  and  z9  and  their  time 

derivatives)  to  functions  of  the  single  variable  9  and  its  time  derivative.  The 
transformation  that  accomplishes  this  is: 
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Zj  -  z^(O)  -  R^# 


(4) 


7. j  =  z9(0)  -  R,?ff 

so  that  the  kinetic  energy  becomes: 

T  =  (^(m^R^2  +  m9R92)  92 

and  the  potential  energy  becomes: 

V  =  mjgz^O)  +  m.2gz.2(0)  +  (m.^-n^R^gfl 

The  Lagrangian  function  is: 


L  =  T  -  V 


and  using  Hamilton's  principle: 


6  f  Ldt  =  0 


We  are  led  to  the  Euler-Lagrange  equation: 

9  = 


liifRi  —  m9R2 


nijRj2  +  in2R92 


In  terms  of  z^  and  z<2,  we  have: 


m2R2  ~  m;Ri 
m2R2  +  mlRl2 


z,  = 


Rl« 


m  R,  -  m2R2 

z9  = - R9g 

rn^Rj  +  m2R92 


(5) 

(6) 


(7) 


(S) 


(9) 


(10) 


(11) 


(12) 


It  is  obvious  that  if  R1  =  R2,  Equations  (11)  and  (12)  reduce  to  Equations  (9)  and  (10)  of 
the  previous  section  [II. l.a]. 

We  assume  that  the  interested  reader  can  integrate  and  interpret  both  the  first  and  second 
integrals  of  these  equations. 


Il  l  c  An  idealization  of  an  Atwood's  machine  wherein  the  wire  connecting  the  two 
masses  is  elastic  is  given  in  Figure  l.c.  The  following  additional  notation  is  used: 
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A  =  cross-sectional  area  of  the  wire 
E  =  Young's  modulus 

f  =  length  of  the  unstressed  wire 

A  P 

k  =  elastic  constant  of  the  wire  = 

Again,  using  the  nomenclature  of  that  figure,  the  kinetic  energy  is: 

T  =  (J)nVi2  +  Q)nV22  (*) 

and  the  potential  energy  is: 

V  =  m1gz1  +  m,2gz2  +  Vs  (2) 

and  where  V  is  the  strain  energy  in  the  wire.  V  is  found,  assuming  Hooke's  law,  as 
follows: 

If  ?  is  the  length  of  the  unstretched  wire,  i.e.  the  length  before  the  masses  and 
m.-,  are  attached,  then  the  strain  after  the  system  is  in  motion  is: 

s  =  (2h  +  ?rk  -  f)  -  (z^+Zg)  (3) 

The  strain  energy  is: 

Vs  =  (l)ks2  (4) 

and  k  is  the  spring  constant  for  the  wire.  Again,  the  Lagrangian  is: 

L  -  T  -  V  (5) 

and  Hamilton's  principle  dictates: 

/l 

6  f  Ldt  =  0  (6) 

h 

The  Euler-Lagrange  equations  lead  to: 

mjZj  +  k(zl  +  z2)  -  k[zj(0)  +  z2(0)]  -  m}g  (7) 

m2z2  +  k(z1  +  z2)  =  k[z1(0)  +  z2(0)]  -  nx^g  (S) 

which  when  integrated  become: 
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Z1  =  m2^  cos  +  B  s’n  <Tt)  +  ~ 


z9  =  cos  crt  4-  B  sin  at)  +  1 


gt2  +  vQt  +  z* 


gt  -v0t  +  z2* 


where  A  and  B  are  constants  of  integration,  Vq  is  an  initial  velocity, 


a2  =  k 


m  1+  m2 
"mlm2 


and 


zl*  +  z2*  ~  7i(°)  +  z2^  " 


'  2  m  1  m2 

r /•  4.  /*i 

‘l+  ‘2 

m  1+  m2 

AE 

and  where: 


(9) 

(10) 


(11) 

(12) 


4*  =  length  of  the  unstretched  wire  as  measured  from  the  top  of  the  pulley  (point  P)  to  the 
point  where  the  mass  nn  is  attached  (i  =  1,2). 


We  infer  from  the  above  equations  that  the  motion  of  the  two  masses  and  m2  in 

an  Atwood's  machine  in  which  the  two  masses  are  connected  by  an  elastic  wire  consists  of 
two  "normal"  modes.  The  first  is  the  motion  that  would  occur  if  the  wire  were  inelastic. 
The  second  may  be  interpreted  as  each  mass  executing  simple  harmonic  motion,  the  two 
harmonic  motions  being  in  phase,  with  the  amplitude  of  the  motion  of  mass  m^  being 

proportional  to  the  mass  m9  and  the  amplitude  of  the  motion  of  the  mass  m2  being 
proportional  to  m^.2 


II.  1  d  With  Friction  on  the  Bearing.  In  this  section,  the  equations  of  motion  for  an 
Atwood's  machine  for  which  there  is  friction  on  the  bearing  are  derived.  Referring  to 
Figure  I  d,  T^  is  the  tension  in  the  wire  that  joins  the  mass  to  the  pulley  and  T2  is  the 

tension  in  the  wire  that  joins  m2  to  the  pulley.  Also  from  the  figure,  it  is  clear  that: 

Z1  +  z2  =  ^  +  ~  *  =  constant  (1) 

where  T  =  length  of  wire  joining  the  two  masses.  Differentiating  Equation  (1)  with  respect 
to  time  yields: 

Zj  +  z9  =  0  (2) 

and 

zi  +z2  =  0  (3) 
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The  force  acting  on  is: 


T1=m1g  +  m1z1  (4) 

and  on  m9: 

T2  =  m2g  +  m2z2  (5) 

If  m„  >  mp  so  that  the  direction  of  motion  is  as  depicted  in  the  figure,  i.e.  clockwise, 
the  forces  balance  out  as: 


+  F 


B 


(6) 


where  Fg  is  the  force  due  to  friction  on  the  bearing.  Denoting  the  coefficient  of  friction  by 
/i,  the  total  frictional  force  is  assumed  to  be: 


Thus, 

or 


FB  “  ^T1  +  T2^ 

(?) 

Equation  (6)  becomes,  using  Equations  (3),  (4),  and  (5), 

T2  =  Ti  +  n{Tl  +  T2) 

(8) 

(mj— n^)  +  /i(m1+m2) 
z2  “  (n^+rr^)  +  ^(m^— m2)  s 

(9) 

If  the  coefficient  of  friction  of  a  particular  bearing  is  unknown,  n  may  be  found  by 
solving  for  it  in  Equation  (9).  A  simple  experiment,  using  known  masses  and  m2,  a 

known  value  of  g  and  determining  z2,  enables  one  to  calculate  /r. 


II.2.a  An  Atwood’s  pendulum  is  defined  as  an  Atwood's  machine  in  which  one  of  two 
masses  is  allowed  to  swing  as  a  pendulum  while  the  other  remains  constrained  to  move  only 
in  the  vertical  direction.  The  pendulum  motion  of  the  one  mass  induces  a  varying  tension 
in  the  connecting  wire;  this,  in  turn,  produces  motion  in  the  second  mass.  It  is  shown  that 
this  motion  can  be  made  periodic  if  the  ratio  of  the  two  masses  and  the  dependency  of  this 
ratio  on  the  initial  conditions  are  chosen  as  to  be  prescribed  in  this  report.  If  this  condition 
is  not  met,  the  m'- '  ion  cor  ists  of  the  superposition  of  two  motions.  The  first  is  motion  in 
a  constant  gravlt"  rnal  field  where  the  effective  "gravity"  is  kg;  the  factor  k  is  determined 
explicitly.  TV  v  cud  is  the  periodic  motion  that  is  the  central  theme  of  this  section  of  the 
report.  During  tne  course  of  the  analysis,  the  fundamental  frequency  of  the  periodic 
motion  is  determh  '  *  It  is  shown  to  be  slightly  higher  than  the  frequency  of  a  pendulum 
of  comparabl  lengr  swinging  in  the  earth's  gravitational  field;  the  factor  is  given 
explicitly.  This  work  is  restricted  to  the  extent  that  small  angle  approximations  are 
introduced  initially  for  trigonometric  functions. 

The  geometry  of  the  configuration  studied  is  given  in  Figure  2.  From  this  geometry, 
assuming  that  the  masses  of  the  pulley  and  the  wire  are  negligible  and  that  the  radius  of 
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the  pulley  may  be  neglected  also,  the  differential  equations  describing  the  motion  of  the 
mass  m«2  are: 


p()  +  2 p9  +  g  sin  0=0 


fl) 


(m^  +  m2)  P  -  m^pO2  +  m^g  -  rn^g  cos  9  =  0 

(2) 

where  dots  indicate  differentiation  with  respect  to  the  time  t.  Equations  (1)  and  (2)  may 
be  simplified  slightly  if  we  observe  that  by  letting  p  =  -ig,  both  equations  contain  g  as  a 
factor  and  we  may  divide  through  by  it.  This  is  equivalent  to  choosing  units  so  that  g  =  1. 
Furthermore,  by  using  the  small  angle  approximation,  these  equations  then  become: 

iff  +2*  0  +0=0  (3) 

(n^  +  m9)i  -  m2<*02  +  (m^  -  m2)  4-  A  =0  (4) 


Case  i)  Constant  * 

We  prove  now  that  there  is  only  one  solution  of  Equations  (3)  and  (4)  for  which 
-*=  -*q  =  constant.  Actually,  we  shall  show  that  assuming  that  t=  constant  leads  to  a 

contradiction.  Using  this  assumption,  Equations  (3)  and  (4)  become 

4Qff+0=  0  (3.1) 

m2V2  ~  (mi  “  m2)  "  \  ~  0  (4.1) 

Differentiate  Equation  (4.1)  with  respect  to  the  time  t  and  divide  through  by  m2  to  get: 


2*Q  M  -  00  =0  (5) 

If  0  =0,  then  Equation  (3.1)  implies  that  ^=0  and  Equation  (4.1)  then  implies  that 
m1  =  m9,  the  classic  Atwood's  machine  case  with  equal  masses.  For  this  classic  example 

there  is  indeed  a  solution  i=  4Q,  a  constant.  If  6  $  0,  then  Equation  (5)  implies  that 

2*0  0=9  (6) 

Equations  (3.1)  and  (6)  are  incompatible.  There  is,  therefore,  no  solution  of  Equations  (3) 
and  (4),  other  than  the  classic  Atwood's  machine,  for  which  *,  is  constant. 

We  proceed  now  to  discover  the  relationship  between  nij  and  m2  for  which  *(t)  and 
9(t)  are  periodic  solutions  of  Equations  (3)  and  (4). 
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Case  ii)  The  General  Case. 

The  perturbation  technique  requires  that  Equations  (3)  and  (4)  be  modified  and 
rewritten  as  follows: 


2  -  (nij  -  m^)  -  -  m,^2  j  (7) 

*$+8  =  -  2*6  (8) 

where  we  have  introduced  a  parameter  €  and  written  =  m^/(mi  +  mg)  and 
no.,  =  m9/(m^  +  mg).  In  Equations  (7)  and  (8),  we  change  the  time  variable  t  to  a 

pseudo-time  t  by  means  of  the  formula  t=  <Ju>  t  and  get 

^'  =  e[mg^0/2  +  w(mg  -  m^)  -  I  mgtu^2  ]  (9) 

*0'  +  v0  =  -2*'0'  (10) 

where  the  '  indicates  differentiation  with  respect  to  t.  Assuming  that  and  are 
analytic  functions  of  6: 

■*=  rQ  +  +  e2r2  +  ... 

0=  00  +  601  +  €2^  +  ... 
c rf  —  u>q  +  T  •" 


ml  —  ^0  +  +  6^2  + 


We  introduce  these  expansions  into  Equations  (9)  and  (10)  and  collect  terms  in  like 
powers  of  6.  Sir>™>  e  is  an  arbitrary  parameter,  each  coefficient  of  the  powers  of  €  must  be 
zero.  This  results  in  setting  up  an  infinite  system  of  pairs  of  second  order  differential 
equations  which  can  be  solved  successively  and  recursively.  The  first  set  of  pairs  is 
obtained  by  setting  €  =  0.  The  two  equations  obtained  are: 

r0-=0  (11) 


r 


+  ‘Vo-  2ro% 


(12) 


The  solution  of  Equation  (11)  is  Tg  =  r^,  a  constant,  since  we  assume  as  an  initial 
condition  that  Iq'(0)  =  0.  The  solution  of  Equation  (12)  then  becomes 


0q  =  0QCos(kt) 


(13) 
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where  0g  is  the  initial  angular  displacement,  k2  =  ^q/tq  and  we  have  assumed  that 
0q'(0)  =  0.  It  will  turn  out  that  o>q  plays  no  role  in  the  final  form  of  the  formulas  for  p(t) 

and  fl(t).  In  terms  of  the  physical  variables,  k2  =  g//>(0). 

If  Equation  (9)  is  differentiated  with  respect  to  the  parameter  e  and  then  e  is  set 
equal  to  zero,  we  obtain: 


V  =  HYoV2-^  +  “W1  V)  (14) 

When  the  values  of  #q  and  9q  ,  as  given  from  Equation  (13)  and  its  derivative,  are 
substituted  into  Equation  (14),  and  remembering  that  Tq  =  Tq,  a  constant,  then,  after  some 
algebraic  manipulation,  Equation  (14)  may  be  written  as: 

V  =  +  “2  +  ~  (15) 

Since  we  are  seeking  a  v eriodic  solution  for  r,  the  secular  term  in  the  general  solution 
of  Equation  (151,  the  term  that  would  give  rise  to  a  quadratic  increase  or  decrease  with 
time  in  r^,  will  oe  eliminated  if  we  set: 

^0  =  “2  +  (16) 

Otherwise,  will  contain  a  term  of  the  form  (~pg  4-  +  (i)m2^g2)t2/2  and,  depending 

on  the  relative  magnitude  chosen  for  /Zq  and  (m.^  4-  (ijm^Q2),  increases  or  decreases 

quadraticallv  with  time.  The  solution  of  Equation  (15),  for  the  initial  conditions 
^(0)  =  rj'fo)  =  0,  is: 


T 


(17) 


When  Equation  (10)  is  differentiated  with  respect  to  e  and  then  £  is  set  equal  to  zero, 
we  obtain: 


ro*I  +  uo0i  ”  ri*o  2ri^6 


(18) 


Since  r^,  0^  and  are  known,  we  substitute  their  respective  values  into  the  right-hand 
side  of  Equation  (18)  to  obtain: 


r0^1  +  w0^1  =  +  + 


[(^)m./02u,0]cos(3kt)} 


(19) 
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We  axe  interested  in  obtaining  the  particular  solution  of  Equation  (19)  that  remains 
bounded  with  increasing  time.  i.e.  we  are  interested  in  excluding  the  secular  term.  That 
part  of  the  particular  solution  of  Equation  (19)  that  arises  f’-cm  the  term 

+  ©“2  V^cos(kt^ 


is  of  the  form  at  sin(kt),  where  a  is  a  constant  depending  on  the  several  constants  in  the 
equation,  namely,  Tq,  Uq,  ,  0q,  m^,  and  k  and  whose  explicit  dependency  on  these  is  not 

important.  What  is  important  is  that  unless  a  =  0,  the  function  t  sin(kt)  oscillates  wildly 
and  it  is  this  oscillation  that  must  be  eliminated.  This  is  done  by  choosing 

=  -(tf)miVu;o  (20) 

With  this  done,  the  particular  solution  of  Equation  (19),  with  initial  conditions 
^(0)  =  ^'(0)  =  0,  is: 


9l  =  ~  cos(3kt)]  (21) 

Note  that  if  is  not  chosen  as  in  Equation  (20),  then  the  solution  for  6^  will 

oscillate  wildly  with  time.  We  interpret  this  by  saying  that  the  mathematical  procedure 
fails.  Thus,  in  order  to  retain  mathematical  viability,  we  must  pick  as  given  in 

Equation  (20). 

If  we  were  to  proceed  no  further,  we  would  have  the  following  approximate  solutions 
to  Equations  (3)  and  (4)  (after  setting  e  =  1): 


*=  r0[i  +  v^cos(2kt)  -  ^ 

(22) 

0=  0Q{cos(kt)  +  (-^)(m.2^Q2)[cos(kt)  -cos(3kt)]} 

(23) 

ml  =  “2  + 

(24) 

“  =  “oi1  -  VI 

(25) 

k2  =  wg/ro 

(26) 

We  proceed  to  get  the  next  higher  order  terms  in  the  expansions  for  4,  0,  m^,  and  u>. 


II 


By  equating  the  coefficients  of  €2  from  both  sides  of  Equation  (9),  we  get 


V  =  '[‘Vi  +  "  “2)1 +  “sMo2  +  2roOT 

+  (2?) 

After  substituting  the  several  quantities  already  found  for  their  respective  values  as  given 
in  the  right-hand  side  of  Equation  (27),  it  may  be  rewritten  as: 

r9*  =  (^-^)m22#Q4a^[4  cos(2kt)  +  9  cos(4kt)]  (28) 

and  where,  in  order  to  preclude  the  introduction  of  a  secular  term  in  ^  we  had  set: 

h  =  OW  <29> 

The  solution  of  Equation  (28),  subject  to  the  initial  conditions  ^(0)  =  r9'(0)  =  0,  is: 

r>2  =  -2-  (m.2^Q2)2rQ[25  -  16  cos(2kt)  -  9  cos(4kt)]  (30) 

1213J 

By  equating  the  coefficients  of  e2  from  both  sides  of  Equation  (10),  we  get: 

r0^'2  UQ^2  trl^l  ^1^1  d" 

^r2^0  +  ^2^0  +  2r2  ^0  ^  =  ®  (31) 

After  substituting  the  values  already  found  for  and  their  derivatives 

where  appropriate,  we  get: 

ro02’  +  =  f-2-]m22i ^05co0[163  cos(3kt)  -  291  cos(5kt)J  (32) 

I2 1 4  J 

and  where,  in  order  to  preclude  the  introduction  of  a  secular  term  in  we  ^ave  se^: 

u2  =  (33) 

The  value  of  flg,  subject  to  the  initial  conditions  ^(0)  =  #2  "W  =  ^oun(^  by  integrating 
Equation  (32)  is: 

6 2  =  m22#Q5[66  cos(kt)  -  163  cos  (3kt)  +  97  cos  (5kt)] 


(34) 


Summarizing  the  results  obtained  thus  far,  we  see  that  after  setting  e  =  1,  we  get  the 
following  approximate  solutions  for  -t(t)  and  0(t): 


*(t)  =  rQ[l  +  (-^)m.2^2(cos  (2kt)  -  1)  + 


w02Pq4(25  -  16  cos  (2kt)  -  9  cos  (4kt)] 


(35) 


0(t)  =  0q{ cos  (kt)  +  (j^-)nv,^2[cos(kt)  -  cos  (3kt)]  + 

Im.:>2i9g4[66  cos  (kt)  -  163  cos  (3kt)  +  97  cos  (5kt)]}, 


_2_ 

.17 


(36) 


where 


ml  =  +  (;)“2V  +  04 

"  =  “oi1  “  (S)m2  V  +  <  m)m2^04l 


(37) 

(38) 


_1  _1  _l 

Since  t  =  (u>)  2  t,  -e(t)  4  't(u)  2  <),  fl(t)  4  0(u  2  t ),  i.e.  in  Formulas  (35)  and  (36) 


replace  t  by  u  *  t.  Finally, 


p{t)  =  g*(w  t) 


and 


9{t)  =  9(u  2  t) 


II. 2. b  Atwood's  Pendululm — non-periodic  motion.  The  configuration,  symbols  and 

references  to  equations  are  as  in  the  previous  section,  II.2.a.  Whereas  in  the  previous 

section  we  used  1  -  92  /2  as  the  small  angle  approximation  for  cos  9,  we  shall  in  this  section 
set  cos  9  =  1.  Equation  (3)  remains  unchanged  but  Equation  (4)  simplifies  to: 

(mj  +  m2)  +  (m^  -  rr^)  =  0  (4.1) 

Having  seen  that  for  the  periodic  solution  to  exist,  we  must  have  m^  >  m2  (to  counteract 
the  centrifugal  force  produced  by  the  swinging  motion  of  the  mass  m2),  it  is  reasonable  to 
expect  that  if  is  less  than  the  critical  value  as  given  (approximately)  by  Equation  (37), 
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then  the  length  -t  would  increase  indefinitely  (or  at  least  to  the  extent  the  physical 
apparatus  would  permit),  whereas  if  m^  is  greater  than  this  critical  value,  then  -t  would 

approach,  and  become,  zero.  These  observations  were  confirmed  numerically,  and  in  fact 
the  computations  were  carried  out  without  assuming  the  small  angle  approximations,  by 

integrating  Equations  (1)  and  (2)  (Mittleman3  and  Zeigler4). 

The  special  case  for  which  =  m0  was  studied  extensively  using  numerical 

integration  in  Mittleman3  and  Zeigler4.  In  5,  Jon  Lee  discussed,  using  asymptotic  methods, 
both  the  short  time  and  the  long  time  behavior  of  this  case  and  provided  a  theoretical  basis 

for  the  numerical  experiments  reported  in  3  and  4. 


II. 3. a  Two  Countered-Balanced  Dumbbells.  The  system  to  be  described  is  depicted 

in  Figure  3.  The  two  masses  and  are  connected  by  a  masT.  "ss  rigid  rod  and 

suspended  by  inextensible  wires  from  the  two  pulleys  p^  and  The  wires  then  run 

horizontally  across  the  top  of  the  apparatus  to  the  pulleys  p^  and  p4  and  then  drop  down 

to  the  two  masses  m^  and  m,  which  are  also  connected  by  a  massless  rigid  rod.  If  the  only 

forces  allowed  to  act  on  the  two  masses  nn  and  m9  are  gravity  and  the  tensions  in  the 

wires  and  there  are  no  initial  conditions  on  these  two  masses  to  produce  motion  out  of  the 
vertical  plane,  this  dumbbell  can  move  only  in  the  vertical  plane.  For  the  dumbbell  made 
up  of  the  two  masses  and  m4,  initial  conditions  can  be  chosen  so  that  the  four  points 

P^,  p4,  m^,  and  m4  need  not  be  coplanar  and  in  fact,  the  dumbbell  m^  -  m4  need  not  lie  in 

the  vertical  plane  containing  p3,  p4_  The  mathematical  description  that  follows  is  based  on 

these  assumptions. 

The  kinetic  energy  of  the  four  masses  is: 


i'i 


Z- 

1 


1=1 


(1) 


The  potential  energy  of  the  four  masses  is: 

4 

v  =  X  migzi  (2) 

i=l 

We  introduce  the  following  vector  notation:  for  i  =  1, 2, 3, 4  (where  appropriate) 

r.  =  (x-,  y.,  z.) 

1  v  i’  ^  P  v 

^i  ~  ^ix’  Pjy’  ^iz^ 
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<1  -  Pi  •  -'ll5 

0  =  r  -  r 

12  rl  2 

«34  =  V74 

g  =  (0,0,  g) 

The  equations  for  the  constraints  are: 

<py  (ri  ~  r2)-(r1  -  r2)  =  ^2*^12 

(?3~  r4)  (r3-r4>  =  ^34-^34 

h  Pi  •  ij*  +  p3  •  V  - 

1  1 

Pi  -  i>l)  •  (?1  -  i>l)lJ  +  l(r3  -  P3)  (r3  -  P3)]2  =  Lj  -  d13 

^  P2  '  72l  +  P4  ■  V  = 

1  1 

[*2  "  P2) '  (i?2  "  M*  +  K*4  "  P4)  (*4  -  P4)]2  =  l2  “  d24 
where  and  L2  are  the  lengths  of  the  two  connecting  wires. 

Hamilton's  Principle  is: 

h  .  . 

6  J  F(t;XpX9,...,Z2,Zg)dt  =  0 

t0 

and  where: 

F  =  T  -  V  +  SAi(t)0i(x1,x2>...,z3,z4) 

and  where  the  A.(t)  are,  for  the  moment,  Lagrangian  parameters.  Carrying  out  the  detail? 
of  the  computation  leads  to  the  following  system  of  equations: 
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rl  -  Pi 


m  r  =  -m  g  +  2A  (t)(r  -  r  )  +  A  (t) - r 

11  1  11  2  3 

— ) 

in2r2  =  — m*2®  —  r2)  d-  j 


m 


3r3  =  -m3g  +  2A2(t)(r3  -  r4)  +  A3(t) 


7*  ^  m  -  p, 

V4  =  -m4g  -  2A2(t)(r3  -  r4)  +  .v4(t) - j — 


to  which  we  must  adjoin  the  four  constraint  equations. 

This  second  order  system  of  differential  equations  is  replaced  by  an  equivalent  system 
of  first  order  equations  so  as  to  facilitate  using  numerical  integration;  in  this  case  we  chose 

Runga-Kutta-Fehlberg6. 


•  — > 

rl  =  V1 


m 


rl  -  Pi 


=  -nijg  +  2A1(t)(r1  -  r2)  +  A3(t) - j 


r2  =  v2 


-4  -» 

r  c 


in2^2  =  ~m2®  -  2Al(t)(rl  “  ?2)  +  A2^ - 


r3  “  v3 


^  ^  ^  Ta  "  Pa 

m3v3  =  -m3g  +  2A2(t)(r3  -  r4)  +  A3(t)  — 


•  *4 

r4  =  v4 


m 


4^4  =  -m4g  -  2A2(t)(r3  -  r4)  +  A4(t) 


and  again  we  must  remember  to  adjoin  the  four  restraint  equations. 

These  equations  are  not  as  yet  in  a  form  suitable  for  numerical  integration  since  we 
do  not  know  how  to  handle  the  four  A-(t).  We  proceed  as  follows.  With  time,  t,  as  the 

independent  variable,  differentiate  each  of  the  four  constraint  equations,  twice.  By 
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adjoining  these  eight  equations  to  the  twenty-four  that  contain  r-  and  v-.  we  can  eliminate 

all  time  derivatives.  The  result  of  these  lengthy  and  tedious  algebraic  manipulations  is  a 
linear  system  of  four  equations  in  the  four  Aj(t).  The  numerical  integration  can  now  begin. 


Pick  initial  conditions  for  r- (0)  and  v .  ( 0 )  consistent  with  the  constraint  equations. 

Although  it  is  possible  to  parameterize  these  four  constraint  equations,  this  is  not 
recommended.  Instead,  use  a  Newton-Raphson  procedure.  With  these  values,  calculate 
the  Aj(0).  Enter  into  the  Runga-Kutta-Fehlberg  routine.  As  the  program  progresses 

through  the  several  steps  of  the  R-K-F  routine,  it  will  be  necessary  to  know  (i.e.  calculate) 

intermediate  values  of  the  T  and  v.;  these  will  be  got  from  the  constraint  equations.  Thus, 

even  for  one  step  in  the  R-K-F  algorithm,  to  get  7. (t^ )  and  w(t,),  it  is  necessary  to  make 

several  intermediate  computations  using  the  constraint  equations.  The  number  of  such 
intermediate  excursions  will  depend  on  which  order  R-K-F  you  chose  to  use.  The  details 
of  the  code  are  not  included  in  this  report;  the  programming  and  running  of  it  was  done  by 
Arnel  Pacia  of  WRDC/FIBG  and  would  be  available  from  that  source. 

While  there  is  literature  on  systems  of  differential  equations  with  algebraic 
constraints,  none  seemed  immediately  applicable  to  this  particular  problem  and  the  method 
described,  therefore,  was  invented  for  this  application. 

Several  computer  runs  were  made  using  this  code  and  aside  from  the  general 
oscillatory  patterns  that  we  had  come  to  expect  there  was  one  surprise.  When  we  picked 
equal  lengths  for  the  two  wires  connecting  m^  -  m^  and  m2  -  m^  and  displaced  the 

m^  -  m^  dumbbell  from  its  rest  position  and  released  it  so  that  it  would  swing,  we 

calculated  the  distance  from  the  center  of  the  m^  -  m^  dumbbell  to  the  center  of  the 

-  m2  dumbbell  measuring  from  the  first  center  over  the  top  of  the  apparatus  and  then 

down  to  the  other  center,  we  found  that  this  distance  changed  and  depended  on  the  angular 
orientations  of  the  two  dumbbells.  It  was  always  less,  albeit  the  amount  was  very  small, 
than  the  common  lengths  of  the  connecting  wires.  While  this  is  of  no  significance  in  this 
problem,  it  portends  difficulties  if  one  were  to  try  to  formulate  a  problem  for  three  masses 
in  a  row,  again  rigidly  connected,  counter-balanced  by  a  similar  arrangement  and 
connected  by  three  wires  of  equal  length.  It  seems  to  this  writer  that  the  problem  is  akin 
to  the  problem  of  supporting  a  rod  at  three  or  more  points;  this  problem  is  known  to  be 
statically  indeterminate  and  elastic  theory  is  required  for  its  solution. 


II. 3. b  Two  Counter-Balanced  Dumbbells  Connected  By  Three  Wires.  From  the 

discussion  presented  in  the  previous  paragraph,  if  three  wires  are  used  to  join  corresponding 
points  on  the  two  dumbbells,  and  one  of  the  dumbbells  is  set  into  motion,  there  could  be  no 
tension  on  the  center  wire.  This  writer  would  infer  then  that  the  problem  of  three 
connecting  wires  will  yield  results  not  different  from  those  obtained  from  the  problem  with 
two  connecting  wires.  The  evidence  is  as  indicated  above  but  a  more  convincing  argument 
may  be  required. 
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II. 4  One  Dumbbell  Counter-Balanced  By  Weights.  Suppose  the  configuration  discussed 
in  II. 3. a  is  modified  by  eliminating  the  constraint  <p^,  that  is,  by  assuming  that  the  masses 

nij  and  m9  are  not  rigidly  connected,  in  fact,  are  not  connected  at  all.  Under  this 

assumption  the  only  motion  possible  for  these  two  masses  is  vertical,  i.e. 

Xf(t)  =  yj(t)  =  x9(t)  =  y9(t)  =  0.  The  equations  of  motion  for  this  situation  are  trivially 

derivable  from  those  given  above  by  paralleling  the  derivation  without  the  constraint. 

We  have  not  carried  out  the  details  and  thus  have  no  numerical  results  to  report. 

As  yet,  another  possibility  would  be  to  initially  have  m^  and  the  rod  connecting 
them  at  rest  and  impart  initial  velocities  to  m ^  and  i  i.r  The  above  equations  certainly 
suffice;  we  have  not,  however,  carried  out  the  details  and  have  no  results  to  report. 

It  would  be  reasonable  now  to  counter-  balance  the  dumbbell  by  three  or  more 
weights  and  study  the  motion.  Obviously  this  modification  circumvents  the  problem  of 
static  indeterminance  and  might  provide  further  insight  into  the  motion.  This  has  not 
been  done. 


II. 5. a  Counter-Balanced  Rigid  Rod  —  With  Two  Wires.  We  start  with  a  rod  of 

length  l,  uniform  density  and  uniform  cross-section  Aq.  The  center  of  mass  of  the  rod 

is  at  the  point  0.  Consider  this  rod  in  an  inertial  coordinate  system  (x,y,z)  such  that  one 

end  of  the  rod  is  at  the  point  r ^  =  (x^.y^Zj)  and  the  other  end  is  at  the  point 

r9  =  (x9,y2,z9).  If  cos  a^,  cos  cr9,  cos  a ^  are  the  direction  cosines  of  the  rod  with  respect 

to  the  inertial  coordinate  s>stem.  then  (cos  Qj,  cos  tv.,,  cos  a^)  =  (r9  -  v^)/l.  The 
coordinates  of  an  arbitrary  point  on  he  rod  is  given  by  the  vector 

r  =  rj  +  (r9  -  r^)<7  0<cr<  1 

In  the  rod,  set  up  a  right-hand  coordinate  system  (x',y',z'),  with  the  x'  direction 
along  the  length  of  the  rod,  the  y'  and  the  z'  directions  Doth  orthogonal  to  x'  and  to  each 
other,  but  otherwise  unspecified.  The  origin  of  this  system  is  taken  at  0.  If  the  linear  mass 
density  of  the  rod  is  given  by  p(x')  >  0  (in  our  case,  p(x')  =  p q),  then  the  mass  m^,  of  the 

rod  is  given  by 


ml  =  / f(x')  =  Pq1 

and  0  =  J p(x' )  x'  dx' 

and  L  =  fp(x')(x' )3dx'  =  (m.  1 2 )  / 1 2 

1  i. 
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jy ,  Let  the  origin  of  the  (x',y',z') 


and  where  the  integration  is  on  the  interval 

.  L  ~ 

system  have  coordinates  (x.  (,.y.  f,.Zjr)  in  the  inertial  system.  Then  an  arbitrary  point  P  in 

space  may  be  referenced  in  both  systems;  the  relation  between  coordinate  values  in  the  two 
systems  being  given  by 


xlc  +  ailx'  +  a^y'  +  a^z' 


13 


y  y,c  -i-  a.jj.x'  +  a.,2y'  +  a.^z' 


7  =  Z1C  a3lX'  +  a32y  +  a332' 

and  where  A  ~  (a-)  is  an  orf  hogonal  matrix.  For  points  on  the  rod.  y '  =  z '  =0.  and  we 
Lave 

x  =  xlc  +  ailx' 

y  =  >'ic  +  a2r 

z=s2lc  +  a3ix' 

and,  in  terms  of  the  angles  previously  defined: 

a,  j  =  cos  Oj 


a0^  =  cos 


a31  =  cos  "3 


To  calculate  the  kinetic  energy  of  the  rod.  we  consider  first  an  element  of  length  dx'. 
The  mass  dm  associated  with  this  element  is  p(x')dx'  —  /?gdx'  and  the  square  of  the 


velocity,  v  .  is: 


and  where 


V-  —  X‘  +  V  +  Z“ 


x  _  xlc  +  allx' 


v  —  V.  +  a..  .X' 

-  •  1  (  j.  i 


z  =  zic  +  a.31x' 
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The  kinetic  energy  of  this  element  is  equal  to  t  l)[n(x' )dxTv;  and  after  summing  over 
the  length  of  the  rod,  the  kinetic  energy  is 


T.  =  hhd 


•  lc 


Te¬ 


rn,  +  !(a 


11 


4-  a21“  +  I] 


This  is  nothing  more  than  the  well  known  result  that  the  total  kinetic  energy  of  a  rigid 
body  is  equal  to  the  translational  kinetic  energy  of  the  body  as  if  all  its  mass  were 
concentrated  at  the  center  of  mass  and  moving  with  fne  velocity  of  the  center  plus  the 
kinetic  energy  of  rotation  about  the  center  of  mass. 

The  three  variables  a^.a^-a.^  are  not  independent  but  are  related  by  the  equation 

a  -,+a,n4a3i3  =  1.  Independent  variables  o  may  i  e  introduced  as  follows:  let 
—  cos  Oy  a,  1  -  cos  dj  sin  Oy  a9,  =  sin  c  ■  sii.  Oy  Then 


•  '■)  y  ,  0  A  l  •  p.  *■  r) 

a,  j -  r  a0j  -r  a  —  0 ^ -  +  sin~p  ~ 

The  potential  energy  of  Tie  rod  is: 

V]  -  f  p(x')  z  gdx'  =  j  p(x')( zJc  f  a.jjX')  gdx' 

=  / /'(*')  dx'  +  a31g  J p<K')x'dx'  =  nijg  z]c 

The  Lagrangian  function  for  tiie  rod  is: 

Since  a  second  rod  is  to  be  included  in  the  configuration  being  studied,  we  would 
repeat  the  above  computation  with  all  subscripts  1  replaced  by  subscripts  2.  This  would 
lead  to  the  following  values  for  T9,  V9  and  L9. 


T2  =  J!’x-n2  +  >2c2  +  V)rn2  +  J^ll2  +  fe2l  ’  +  h3l  >  ‘2 


and  the  direction  cosines  la..)  have  been  replaced  by  the  corresponding  direction  cosines  for 
the  second  rod  (im).  The  potential  energy  is  given  by: 


and  the  Lagrangian  fun.*  lion  -  IT  V9.  There  remain  the  constraints  imposed  by  the 
wires  connecting  the  two  rods. 

We  d'welop  Tie  fa  r  for  two  vires  connecting  the  two  rods.  Refer  to  Figure  4.  The 
positions  of  the  four  puilqw  in  the  inertial  system  are  given  by  Pj  =  (p^,p9j,p3j), 
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i  =  1 .2,3.4.  We  assume  that  all  four  p.^  are  at  the  same  height  above  tha  ground  plane, 
although  this  is  not  essential.  The  distances  !pj  "  P3!  =  |p9  ~p4|  =  d^  -  d94  and 
( pL  -  p0:  =  ( P3  -  P4i  =  dpV  =  d,^4  are  assumed  to  be  constant.  We  shall  refer  to  the  wire 


passing  over  pulleys  1  and  3  as  wire  1  and  the  v.  ire  passing  over  pulleys  2  and  4  as  2  and 
denote  the  length  of  wire  1  by  and  the  length  of  wire  2  by  ly.  Wire  1  is  connected  to  rod 

1  at  the  point  (x,*  ,0.0)  in  the  coordinate  system  fixed  in  the  rod  1  and  is  connected  to  rod 

2  at  the  point  (x  *  .0.0)  in  the  coordinate  system  fixed  in  that  rod.  In  the  common  inertial 
coordinate  system,  the  coordinates  of  these  two  points  are: 

'•xlc  +  ailxf  ’ylc  +  a2ixr  Tc  +  a31xl*  ■  and  -x2c  +  bllx3*  ,y2c  +  b21x3  ,z2c  +  b31' 
res  pec  lively.  We  shall  refer  so  this  first  triplet  of  values  by  a,  and  the  second  triplet  by 


*  \ 

3  > 


'W 


Tiie  length  /  of  wire  1  is  then  given  by 


Ff.  I QL  -  Py  {  +  d13  +  iP3  -q3l  =  \ 

In  strictly  analogous  fashion,  wire  2  is  connected  to  rod  1  at  the  point  (x^*  ,0,0)  and  to  rod 

2  at  the  point  (x£  .0.0).  In  the  inertial  coordinate  system,  we  would  have 

P9  =  bxlc  al  lx2  ’  ylc  "h  a,21x2  '  zlc  "h  a3lx2  ^  and 

q4  =  (x2c  +  b  1 1  x4*  '  y2c  +  b2ix4*  ’  z2c  +  b.3 1 x4  )•  Tbe  len°th  l2  of  wire  2  is  then  §iven 

by:  ~  “  ‘  . 


iq.2-p2l  -r  d24  +  j  P4  -  q4  |  =  l2 


Note  that  T ^  and  T9  represent  the  constraints  on  the  system.  The  coordinates  of  the  pj  are 
assumed  to  be  fixed  but  the  coordinates  of  the  cjj  are  variable. 

While  the  time  t  is  the  only  independent  variable,  the  dependent  variables  are  x^c, 
v^r,  Zjr.  9 j,  d». ,  x0  ,  y9(;,  z9c,  #9,  60-  We  are  now  in  a  position  to  form  the  generalized 
I.agrangian: 

L*  —  I, j,  +  L9  +  A|(t)  +  A9(t)  r9 

and  obtain  the  Euler-Lagrange  equations  by  taking 

6  J  L*  dt  =  0 
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We  leave  the  computation  at  this  point  because  it  parallels  ix<:  n  ’  . dually  and 

numerically  the  work  presented  above.  Computer  codes  were  not  written  for  this  case; 
time  did  not  permit. 

II. 5. b  Counter-Balanced  Rigid  Rods  —  With  Mojo  Than  Tv:>  Wires.  If  the  two 

wires  are  attached  symmetrically  to  the  rods  at  tire  points  AO),  and  (±>:^  .0,0)  and 

rod  2  is  set  swinging  while  the  initial  conditions  for  rod  l  are  chosen  so  that  it  can  move 
only  in  the  [yrz'.j  plane,  then.  I  believe  that  the  distance  between  the  centers  of  the  two 

rods,  as  measured  over  the  top  of  the  apparatus,  would  be  less  than  the  distance  similarly 
defined  if  measured  when  both  rods  are  at  rest .  Thus,  for  example,  if  a  third  wire  looped 
over  the  top  of  the  apparatus  were  to  join  ’he  two  c r  .iters,  umrc-  would  be  slack  during 
most  of  the  time  that  rod  2  is  swinging.  The  situation  reeuces  then  to  one  of  static 
indeterminancy  and  it  is  questionable,  in  my  mind  that  Cm  method  of  approach,  namely 
two  counter-balanced  rigid  rods  connected  by  three  or  more  v,  i-  '  can  be  correctly  solved 
by  simply  adding  additional  constraints,  jo--.  would  need  be  m<..  :>•  to  elastic  theory. 

If  one  forgoes  the  second  rod  and  attaches  t  he  one  rod  to  a  set  of  isolated  point 
masses  by  wires  stretched  over  the  top  of  the  apparatus  the  problem  becomes  tractable 
and  the  methods  described  above  can  be  used  to  analyze  and  describe  the  motion.  These 
details  have  not  been  carried  out. 


II. 6. a  Two  Counter-Balanced  Flexible  Rods  —  With  Two  Wires.  In  this  section,  we 

set  up  the  Lagrangian  function  of  two  counter-balanced  flexible  rods.  The  general 
appearance  of  the  apparatus  is  as  for  the  case  of  two  rigid  rods;  as  may  be  expected, 
however,  further  analysis  is  required. 

We  start  with  an  inertial  coordinate  system  (x^.y^.z^)  and  consider  a  rod  of  uniform 

density  and  cross-section  whose  center  of  mass  is  at  the  point  0;  (x^y^.Zj0)  in  the 

inertial  coordinate  system.  Let  y9,zvj  be  a  coordinate  system  whose  origin  is  at  0  and 

whose  axes  are  and  will  remain  parallel  to  those  of  the  (x,  ,ypZ^)  system.  The 

axes  will  move  as  the  center  of  mass  of  the  rod  moves.  Let  (x3,y3,z3)  be  a  coordinate 

system  whose  origin  is  also  a  0  but  whose  axes  coincide  with  the  principal  axes  of  inertia  of 
tire  rod;  the  axis  lying  along  the  long  axis  of  the  rod.  (Figure  5. a.) 


Let  P  be  a  point  m  mace  whose  coordinates  relative  to  the  (x.^y.^.z^)  axes  are 
and  whose  coordinates  in  the  (x.,.y9  z.vj  system  are  (o-.v. then 


(o2.T2.72)  =  A(aij)(mrB3!73)1 

where  A(a-)  is  an  oithogona!  matrix  specifying  the  orientation  of  the  (x.^.y^.z.^)  system  in 
the  (x9.y,vz9)  ~yst cm  and  ( A3A3,73)r  stands  for  the  transpose  T  («3./?3,73). 


If  the  coordinates  of  the  point  P  in  the  inertial  system  are  (a^3^, 7^),  then 

=  (xi°>y  i°’zi°)  +  A(aijj(«3^3,73)T 

We  consider  now  the  rod  flexed  about  the  >3  direction  Before  bending,  a 
cross-section  in  a  plane  perpendicular  to  the  y3  direction  remains,  after  bending,  in  a  plane 
perpendicular  to  the  >3  direction.  Such  a  cross-section,  in  a  plane  parallel  to  the  (-'<3123) 
plane,  is  depicted  in  Figure  o.b. 

Suppose  now  that  the  rod,  before  being  bent,  had  been  sectioned  into  small  lengths 
by  a  set  of  planes  perpendicular  to  the  x3  axis.  Let  us  look  at  one  such  section  and  the 

partitioning  planes  after  bending.  Grossly,  this  is  indicated  in  Figure  5  c;  the  two  planes, 
which  had  been  parallel  before  bending,  now  intersect  in  a  line  parallel  to  the  y3  axis  and 

pass  through  the  point  marked  0"  in  that  fgure.  An  enlarged  version  of  that  section  is 
given  in  Figure  5.d. 

If  (03, $3, 73)  are  the  coordinates  of  the  center  of  mass  of  the  element  of  volume 

depicted  in  Figure  5.d,  then  we  set  up  a  fourth  coordinate  system  (x,,y  ,,z.)  whose  origin  is 

,y3,z3)  axes.  In  this  coordinate  system, 

the  mass  and  volume  of  the  section,  which  we  denote  by  Am  and  Av,  are  related  by  the 
equation  Am  =  PgAv. 

Let  4,74)  be  the  coordinates  of  a  point  P*  in  the  element  Av  depicted  in  Figure 
5.d.  Since  the  (x4,y4,z4)  system  is  parallel  to  the  (x3,y-3Z3)  system  and  (0-3, $3, 73)  are  the 
coordinates  of  the  center  of  mass  of  the  element  Av  in  the  system,  then  the 

coordinates  of  P*  in  the  inertial  system  are: 

(x^y^)  =  (xjVj0^0)  +  A(ajj)(Q'3+a4,/?3+/?4,73+74)T 
The  kinetic  energy  of  this  element  of  mass  is: 

1  // /V’T+V+h2*  dx4dy4dz4 

(Av) 

where 


at  (03,^3,73)  and  whose  axes  are  parallel  to  the  (x3 


(Xpy^Zj)  =  (x10,y10,z10)+A(ajj)(a3,^3,73)T+A(ai.)(a3+a4,/?3+/?4,73+  74)T 
The  "phase  space"  variables  are: 


°,z1°),  A(ajj),  (or3,/?3, 73).  («3^3,73) 
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since  the  variables  are  dummy  variables  in  the  integ!a,;.„i 

To  calculate  the  potential  energy  of  this  element  of  mass,  v,e  need: 

Zl  =  Zf9  +  a31(a;3-fo1)  +  a32(  J;}+^4)  +  a3  ,  .  F  y}) 

so  that  the  potential  energy  is: 

/ / /^Ogfzl°  +  a31(o3  +  r>4)  4  a32^3+J4'  +  a33( 73'r  -1-  dx4d>4dz4 
(Av)  . ' 

=  P0§lzq°  +  a3 1  ^ °'3 )  ^39(^9!  T  a33-  73)]  <A'r 

=  ^.gJZ]0  +  a3l(o-3)  +  a33(;?3)  +  ^33t  73)] 


The  "phase  space"  variables  arc.  z^0,  t3),  and  (a31,a39,a33). 

We  calc’daie  the  strain  energy  U  in  the  rod  as  follows:  Let  (dU/dv),  (where  dv  is  an 

element  of  volume),  be  the  work  of  deformation  so  that  (dU/dv)  =  1  ere  where  a  is  the 

stress  in  the  x3(=x4)  direction  and  6  the  corresponding  strain  (deformation).  From 

Hooke's  law  o  -  Eg.  where  E  is  Young's  modulus,  and  so  the  quantity  we  need  to  calculate 
2 

is  I  ^7-  =  A  Ee2.  We  refer  again  to  Figure  o.d. 

Prior  to  bending,  the  sections  mm  and  pp  were  parallel;  after  bending  these  two 
sections  intersect  in  a  line  (whose  trace  in  Figure  o.d  is  O')  parallel  to  the  y3  axis.  The 

longitudinal  fibres,  (those  measured  along  the  x^  axis)  undergo  extension  on  the  convex 

side  of  the  neutral  axis,  whereas  those  on  the  concave  side  are  compressed.  The  arc  nq  is 
the  trace  of  the  surface  in  which  the  fibres  are  neither  compressed  nor  extended  during 
bending  (the  neutral  surface). 

The  elongation  of  a  fiber  at  a  distance  above  the  neutral  surface  may  be  found  by 

first  drawing  qs  parallel  to  mn.  We  look  next  at  the  two  triangles  stq  and  nqO'.  The 
radial  line  O'qt  r  uts  the  two  circular  arcs  at  right  angles,  i.e.  we  assume  that  the  angles 
nqO'  and  stq  are  both  right  angles.  Also,  since  sq  is  parallel  to  nO' .  the  line  tqpO'  may  be 
considered  a  transversal  intersecting  two  parallel  lines  and  thr. '  igle  sqt  =  angle  nO'q. 
The  two  triangles  are  similar  and  their  corresponding  sides  are  proportional.  The  "unit" 
extension  of  the  fiber  rs  is  (st/nq)  =  (qs/O'n)  =  ~(z  Jp),  where  p  =  radius  of  curvature  of 

the  arc  nq.  (Note:  do  not  confuse  p  and  p^,  the  density  of  the  rod;  the  negative  sign  comes 

from  the  fact  that  in  elementary  calculus,  the  curvature  is  considered  positive  if  the  curve 
bends  convex  upward.)  Thus,  the  strains  of  the  longitudinal  fibers  arc  proportional  to  their 
distances  from  the  neutral  surface  and  inversely  proportional  to  the  radius  of  curvature  of 
the  neutral  surface.  Since  we  are  assuming  that  (dz3/dx3)  <<  1,  we  are  assuming 


that  the  curvature  may  be  taken  as  (d2z3/dx32).  This  permits  us  to  write  the  unit  energy 


(dU/dv)  =  lE[z4(d2z3/dx32)]2 


Thus: 


(AU)x  =  iE(d2z3/dx32)2  f  f  f  z42dx4dy4dz4 


lE(d2z,/d::2)2(A::  )  f  f  z,2dv.dz, 
2  V  S'  O  '  v  4 '  J  J  4  '  •*  -i 


If  we  let 

1  =  //  Vd-v4dz4  ’ 

then 

(AU)x=iE(d2z3/dx32)2(Ax4)I 

=  iEI(d2z3/dx32)2(Ax3) 

Continuing,  the  kinetic  energy  of  the  rod  = 

/ tyofff  (Xi2+y12+zi2)dx4dy4dz4]  dx3  , 

(Av) 


the  potential  energy  of  the  rod  = 

/  Amg  [Zl°  +  a31(o3)  +  *32{J3)  +  a33(73)]  dx3 
and  the  strain  energy  of  the  rod  = 


/  ^EI(d2z3/dx32)2(dx3). 

The  integration  in  the  x3  direction  is  over  the  length  of  the  rod,  {-~  ,  i]. 

This  computation  must  be  repeated  for  the  second  rod  which  is  to  be  used  to 
counter-balance  the  first  one.  In  addition,  and  paralleling  the  method  used  in  Section  II. 5, 
two  constraints  and  T^  are  to  be  introduced.  The  Lagrangian  function  takes  the  form: 

L  =  kinetic  energy  -  potential  energy  -  strain  energy  +  A4(t)T4  +  ^(t)!^ 

The  Euler-Lagrange  equations  are  then  obtained  from  Hamilton's  principle: 
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L  dt  =  0 


c0 


Whereas  in  the  previous  examples,  the  Euler-Lagrange  equations  were  ordinary 
differential  equations  with  algebraic  constraints,  in  this  case  they  are  partial  differential 
equations  with  constraints.  Numerical  procedures,  such  as  the  finite  element  method, 
would  have  to  be  extended  to  take  into  account  these  constraints. 


II. 6. b  If  more  than  two  wires  are  used  to  join  the  two  rods,  the  problem  assumes 

another  dimension  of  complexity.  In  the  previous  section,  the  flexibility  was  assumed 
limited  to  vibrations  about  the  y^  axis;  the  bar  was  treated  as  if  it  were  rigid  with 

respect  to  bending  in  the  direction.  Whether  a  third  wire  would  provide  support  or 

simply  be  slack  when  the  rod  is  vibrating  only  about  the  y^  direction  is  not  known. 

In  order  to  get  a  better  understanding  of  what  is  happening,  it  may  be  necessary  to 
remove  the  restriction  that  the  connecting  wires  be  inextensible  and  allow  for  elastic 
connecting  wires  as  was  done  in  Section  II.l.c. 

The  comments  of  this  section  are  purely  speculative.  It  may  be  possible  to  get 
theoretical  answers.  Numerical  results  should  be  obtainable  with  reasonable  effort. 


26 


References 


[1]  Goldstein,  Herbert.  Classical  Mechanics.  Addison-Wesley  Press,  Inc.,  Cambridge, 
Mass.,  1950,  p.  15. 

[2]  Mittleman.  Don.  The  Elastic  Atwood's  Machine.  AFWAL-TM-S6-197-FIBG, 
Wright-Patterson  AFB,  OH,  45433. 

[3]  Mittleman,  Don.  Large  Space  Structures  Dynamic  Testing,  a  report  submitted  to: 
Southeastern  Center  for  Electrical  Engineering  Education  under  contract  to  the  Air 
Force  Office  of  Scientific  Research.  Flight  Dynamics  Laboratory,  Wright-Patterson 
AFB,  OH,  45433,  20  July  1984. 

[4]  Zeigler,  Michael  L.  Low-Restraint  Suspension  Space  Structure  Dynamic  Testing, 
AFWAL-TM-S5-196-FIBG,  Flight  Dynamics  Laboratory,  Wright-Patterson  AFB, 
OH,  45433,  April  1985. 

[51  Lee,  Jon.  Counter-Balanced  Pendulum.  AFWAL-TR-S6-3014,  Wright-Patterson 
AFB,  OH,  45433. 

[6]  Burden,  R.L.,  Faires,  J.D.,  Reynolds,  A.C.  Numerical  Analysis.  2nd  Edition, 
Prindle,  Weber  k  Schmidt,  Boston,  Mass.,  p.  216. 


27 


Figure  l.a 

The  classic  Atwood’s  machine. 


Figure  s.c 

A  section  of  a  slice  of  the  rod, 
a  cross-section  perpendicular  to  the  y3  axis, 

after  bending. 
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